#### trained on right trolls, right trolls only
pdf("figs/FigureA3_maga_politicization_replicate_psrm.pdf", width=4.1, height=5)
par(mfrow=c(1,1), mar=c(6,3,4,2), mgp=c(4,1,0))
reference <- na.omit(with(
    subset(
        aggregated_tweets,
        tweet_date >= "2016-01-01"
    ),c(sorting1rr)))
with(
    subset(
        aggregated_tweets,
        tweet_date >= "2016-01-01"
    ),
    plot(
               y=rep(as.Date(tweet_date), 1), x=c(sorting1rr), type="n",
        col="red",
        bty="n",
        ylab="Tweet Month",
        xlab="\n Dimension 1 (scale in sd's)",
        lwd=3, cex.lab=1.3,
        main="Right trolls only",
        xlim=c(-1.5, 1.5),
        las=1,
        yaxt="n"
    )
)
axis(
    2,
    seq(as.Date("2016-01-01"),as.Date("2016-11-01"), by="quarter"),
    labels=gsub(
        " ", "\n",
        as.yearmon(
            seq(as.Date("2016-01-01"),as.Date("2016-11-01"), by="quarter")
        )
    ),
    las=2
)
with(
    subset(
        aggregated_tweets,
        tweet_date >= "2016-01-01"
    ),
    lines(
        y=as.Date(tweet_date), x=sorting1rr,
    type="l",
        col="red",
        lwd=5
    )
)
abline(v=0, lwd=3, lty=2)
abline(v=c(-0.5, 0.5), lwd=1, lty=3)
with(
    subset(
        aggregated_tweets,
        tweet_date >= "2016-01-01"
    ),
    mtext(
        text=round(max(reference), 1), side=4,
        at=as.numeric(as.Date(tweet_date)[which.max(sorting1rr)]),
        las=1, adj=-0.5, cex=0.8
    )
    )
## ##
dev.off()
